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The phase space F of quantum mechanics can be viewed as the complex 
QQ , projective space CP" endowed with a Kahlerian structure given by the Fubini- 

\ Study metric and an associated symplectic form. We can then interpret the 

Schrodinger equation as generating a Hamiltonian dynamics on F. Based 



O 



upon the geometric structure of the quantum phase space we introduce the 
corresponding natural micro canonical and canonical ensembles. The resulting 
density matrix for the canonical T-ensemble differs from density matrix of the 
conventional approach. As an illustration, the results are applied to the case 
CN ' of a spin one-half particle in a heat bath with an applied magnetic field. 

OO . PACS Numbers : 05.30.Ch, 02.10.Rn 
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^ ; I. QUANTUM PHASE SPACE 

^ ' With a view to certain geometrical generalisations we reformulate the conventional quan- 

tum theory in the following manner. Consider a finite dimensional complex Hilbert space 
^2 ■ T-Cq'^ with complex elements Z"". We identify the state vector Z°' with its multiples AZ" 
^ ! (A e C — {0}) to obtain the complex projective space CP". Here we use Greek indices 

qh| (a = 0, 1, ■ ■ ■ , n) for vectors in 'HJq^, and regard CP" as the true 'state space' of quantum 

mechanics. The complex conjugate of is written Z^, with the Hermitian inner product 
Z°'Zcf We let the state vector represent homogeneous coordinates for CP". Then, the 
complex conjugate of a point P" in CP" corresponds to the hyperplane PaZ"" = 0. The 
points on this plane are the states orthogonal to the original state P". Thus, CP" is equipped 
with a complex conjugation operation that maps points to hyperplanes of codimension one, 
and vice-versa. 

Distinct points X°' and Y" are joined by a complex projective line L"^ = X'^F^', 
representing complex superpositions of the original two states. The quantum mechan- 
ical transition probability between two states X" and Y" is the cross ratio k = 
X^YaYl^ Xj^l X'' XryY^Y^. If the system is in the state Y" and a measurement is made 
with the projection operator X°'Xp/ X"* X^., then k is the probability that the result is the 
eigenvalue unity. The angle 6 defined by /? = cos^{9/2) can be interpreted as the distance 
between the states X" and F". Then if we set 9 = ds, X° = Z", and = dZ"", 
retaining terms to second order, we obtain the unitarily invariant Fubini-Study metric 
on the state space, given by 

ds^ = 8(Z"Z„)-2z[°dZ^lZ[«rfZ^] . (1) 
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Now suppose we write for the Hamiltonian operator, assumed Hermitian. Then, if we set 
h = 1, the Schrodinger equation reads dZ°' — iH^Z^dt. However, the projective Schrodinger 
equation 

Z^'^dZ^^ = iZ^'^H^^Z'^dt , (2) 

which ehminates the superfluous freedom associated with the scale of Z", is defined directly 
on the state space. 

An alternative way of viewing this structure is to regard CP" as a real manifold F of 
dimension 2n, equipped with a Riemannian metric Qab and a compatible symplectic structure 
VLabi satisfying g"'^Q.ac^bd = Qcd and Va^^tc = 0. We use Roman indices (a = 1, 2, ■ ■ • , 2n) 
for tensorial operations in the tangent space of F. The Schrodinger equation then takes 
the form of a Hamiltonian flow — 2QF'^VbH on F. Here the generating function H{x) 
is given at each point x in F by the expectation i7^n^(a;) of the Hamiltonian operator, 
where n^(a;) = Z^ZajZ^Z^ is the projection operator for the point x E F represented by 
the state vector Thus, in quantum mechanics the dynamical trajectories are given 

by a symplectic flow on the quantum phase space F, generated by the Hamiltonian function 
H{x). The flow is necessarily a KiUing fleld, satisfying 

V(„6) - , (3) 

where = 9abC^- In other words, the isometries of the Fubini-Study metric on F can be 
lifted to TC^^ to yield unitary transformations. Conversely, given the Killing flcld we recover 
the observable function H{x) on F, up to an additive constant, by use of the relation 

^^"'V„6 = 2in+l)iH-H) , (4) 
where H = H^/{n + 1) is the average of the eigenvalues of Hg. 



II. QUANTUM MICROCANONICAL POSTULATES 

In this paper we use the metrical symplectic geometry of the quantum phase space as 
the basis for a reexamination of the traditional hypotheses of quantum statistical mechanics. 
To this end, we observe that the quantum phase space F admits a natural foliation by 
the {2n — l)-dimensional hypersurfaces £e determined by level values H{x) = E of the 
Hamiltonian function. The number of quantum mechanical microscopic conflgurations (pure 
states) with expected energy in the small range £^ to £■ + AE is then given by 0,{E)AE, 
where the state density Q{E) for energy E is 

Here da"- = g°'''ebc---ddx'^ ■ ■ ■ dx'^ is the natural vector- valued (2n — l)-form on F. In the case of 
an isolated quantum mechanical system with energy in the given small range, we can adopt 
the notion of the microcanonical ensemble, and identify the entropy by use of the Boltzmann 
relation S{E) = \n{Q{E)AE). Here, we implicitly assume what might be called the quantum 
microcanonical postulate, which asserts that for an isolated system in equilibrium all states on 
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a given energy surface in the quantum phase space are equally probable. As a consequence, 
the temperature of such a system is given by /? = dS{E)/dE, where f3 = 1/kT. Thus for 
an isolated quantum system with expected energy E the equilibrium configuration is given 
by a uniform distribution on the energy surface £e, with entropy S{E) and temperature 
T{E) as given above. The corresponding probability density function on F, which we call 
the microcanonical T-distribution, is 

= 5{H{x) - E)/n{E) . (6) 

It is a straightforward exercise to verify that the state density Q{E) = /p(5(if(a;) — E)dV, 
which appears here as a normalisation factor, is consistent with expression (|^). 

A general measurable function F{x) on F represents a nonlinear observable in the sense 
of |jl|J^. We interpret F{x) as the conditional expectation (F)^. of the observable F in 
the pure state x. The unconditional expectation of F in the microcanonical T-ensemble is 
then given by {F)e = Jp F{x)^E{x)dV. In the case of a conventional linear observable we 
have F{x) = F^n^(x). Then the unconditional expectation in the state fiE^x) is {F)e = 
F^^'p{E), where the quantum microcanonical density matrix fi'^{E), parameterised by E, is 

= J^U'^^{x)^iE{x)dV . (7) 

Providing we only consider linear observables, the state is fully characterised by ^'^{E). Now 
suppose W{E) denotes the total phase space volume for states such that H{x) < E. Then 
the density matrix ^i'p{E) can be calculated explicitly by use of the formula 



III. CANONICAL r-ENSEMBLE 

Next, we motivate a construction for the corresponding canonical ensemble, representing 
the situation where a quantum mechanical system is in contact with a heat bath at a fixed 
temperature. The results we obtain are related to those that follow from the conventional 
density matrix for the canonical state, though differ in certain essential respects. First we 
note that the projection operator IV^ijj), for a given state y, can itself be thought of as 
the density matrix corresponding to a mass distribution 5y{x) on T, concentrated at the 
point y. Therefore, we can write I^iy) = JpIl'^{x)Sy{x)dV . Then the conventional density 
matrix of quantum statistical mechanics, which we denote r^, is determined by a point mass 
distribution with density 

r{x) = Z{P)-'Y.^-^^%.i^) (9) 

k 

on r, concentrated at the energy eigenstates yk, with Ek being the corresponding eigenvalues 
and Z{i3) = J2k e~^^'' the associated partition function. Hence, for the conventional density 
matrix we have 
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r1 = Z(/3)-^^e-^^'=n^(y,). (10) 

k 

It should be evident that, despite its general usefulness and wide acceptance as a funda- 
mental basis for quantum statistical mechanics, the conventional density matrix is in 
some respects unnatural in the absence of some mechanism forcing the system to energy 
eigenstates, since its phase space distribution is characterised by 5y^{x). Thus we are led 
to consider an alternative expression, suggested by the geometry of T, in which we do not 
assume a concentration on eigenstates. 

Consider two systems I and II and let them make contact to form a combined system 
I+II. We suppose that the systems only interact weakly, so the support of the resulting state 
on the quantum phase space ri+2 of the combined system is concentrated on the product 
A-2 C ri_|_2 of the phase spaces Ti and r2 for the systems I and II. We have in mind the 
situation where system I represents a heat bath with a given inverse temperature /5, while II 
represents a small system immersed in the bath. We wish to calculate the probability density 
p{E2) that system II lies on the energy surface iii -^2, conditional on the total energy of 
the combined system I+II lying in the small range E to E + AE. Now, according to the 
quantum microcanonical hypothesis, the equilibrium distribution over the energy surface £e 
in ri.2 is uniform. Therefore, conditional on a given value E for the combined system I+II, 
the joint probability density function pe{Ei, E2) for system I to lie on Se^ in A and for 
system II to lie on Se2 in A is 

^(FF) KE-E,-E,)n,{E,)n,{E,) 

Here, Q,i{Ei) and Vt2{E2) are the state densities for the energy surfaces Ee^ in A and £e2 
in r2, and the state density Q,i.2{E) for the energy surface £e in r-1.2 is 

/oo 
Oi(£; -e) 02(e) fie . (12) 
-00 

Our goal is to calculate the conditional probability density for system II to have energy £'2- 
This is given by ^(£'2) = J^^pe{Ei, E2)dEi, which imphes 

Since system I, the heat bath, is in equilibrium, by the microcanonical hypothesis it has en- 
tropy Si{Ei) = ln(fii(£'i)A£), and hence temperature P{Ei) = dSi/dEi. If E2 < E, then 
to first order in E2 we have Vli[E—E2) = Qi{E)e^'^^''' , where P{Ei) is the inverse temperature 
of the heat bath. As a consequence we deduce that the probability density ^(£2) for system 
II to he on the energy surface in £2 is given by ^(£2) = n2{E2)e-^^^ / J^^ fi2(e)e-^'de. 
More succinctly, now we drop the subscript 2, and regard the conditioning as implicit in the 
specification of the parameter /3. Then for the probability density of the energy distribution 
on r we have 

0(e) exp(-/?e) 

Pie) - , (14) 
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where ^(e) is the state density on F per unit of energy and Z{P) = Q{e) exp{—Pe)de is 
the partition function. By the quantum canonical postulate we mean the assumption that a 
small quantum system in equilibrium with a heat bath will be in a state characterised by a 
distribution over the energy surfaces of F with density (Q), having a uniform distribution 
on each such surface. The corresponding probability density function p/3(x) on F is 

exp{-l3H{x)) 

P,i^) = ^(^) ' (15) 

which we call the canonical F -distribution. Alternatively, we can write P/3(x) = 
J^p{e)fie{x)de, expressing Pjs^x) as a weighted average of microcanonical ensembles, where 
p(e) is given as in The associated density matrix for the canonical T-distribution is 

therefore 

p^(T) = |^n^(x)p^(x)cir . (16) 



IV. TWO-STATE SYSTEM 

To illustrate how the consequences of the canonical T-distribution ([T5|) differ from those 
of the conventional treatment in terms of the density matrix r^, we study the example of a 
spin one-half particle in a heat bath. First we review briefly the geometry of the Schrodinger 
dynamics for a spin one-half particle. The quantum phase space in this case is CP^, which 
can be viewed as a sphere S"^. The poles of the sphere correspond to the energy eigenstates, 
with eigenvalues Eq = —h and Ei = h, and the Schrodinger trajectories are generated by a 
rigid rotation of the sphere about the axis through the poles, the two stationary points. This 
rotation gives rise to a Killing field on the sphere, where the angular velocity is 2h, the energy 
difference between the poles. Thus the Schrodinger evolution generates a latitudinal circle. 
Now we introduce a complex Hilbert space with coordinates Z°' {a = 0, 1) which we regard 
as homogeneous coordinates for CP^. The complex conjugate of is the 'plane' Za in CP^, 
which in this dimension is a point. The point corresponding to Za is = e'^^Zjs, where 
is the natural symplectic form on the two-dimensional Hilbert space. The formalism in this 
case is equivalent to the algebra of two-component spinors. By use of the spinor identity 
2x[«F/31 = e^^X^F^, where X'^eap = Xp, we obtain ds'^ = AZ^dZ^^Z^dZ^ / {Z^Z^f for the 
Fubini-Study metric, and Z^dZ^ = iHapZ'^Z^dt for the projective Schrodinger equation. 
The Hamiltonian has the representation 

= 2h^i^ , (17) 
^ P^py ' ^ ^ 

where and are the stationary points, satisfying H^P^ = hP°' and H^P^ = —hP"'. 
It follows that the velocity of the trajectory through F is 

^ = 2/1 sin ^ , (18) 

a special case of the Anandan-Aharonov relation 0]. Here, 6 is the distance from to 
P°, given by the angular coordinate measured down from the north pole. The transition 
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probability from to the north pole P" is (l + cos^)/2, and for the evolutionary trajectory 
we obtain 

Z° = cos(e/2)e*('^*+'^)p" + sin(e/2)e-*('^*+'^)p" , (19) 

where 6 and are the initial coordinates on the sphere for the state at t = 0. The vector 
Z" is normalised by setting P°'Pa = —PaP"" = 1- The expectation E of the energy for this 
state is E = hcos9, and its variance is h^siri^ 6. 

With this in mind we now examine the situation where the spin one-half particle is 
immersed in a heat bath. The conventional density matrix r^(T) for the system in this case 
can be expressed as a superposition of projection operators onto energy eigenstates, with 
Boltzmann weights: 



' r- 



(20) 



/3 e-/3So _^ g-/3£;i 

This can also be written = Z~^ exp{—i3H^). The fact that has trace unity follows 
from the normalisation condition. For the expected energy H^r^ we then obtain E = 
—htaxih{ph). 

Next we apply the canonical T-distribution to this problem. For this ensemble the 
probability distribution of the energy in the spin one-half case is 

p{E)dE = -27iZ{(3)-^smeexp{-(3hcose)de , (21) 

where E = h cos 6. The normalisation condition on the density function p{E) implies Z{f3) = 
47r(/5/i)~^ sinh(/5/i), from which the expected energy E = —d\nZ/d(3 can be deduced. The 
result is E = kT — fiB coth {fiB/kT). Here we have reinstated Boltzmann's constant k, the 
particle's magnetic moment fi, and the external magnetic field strength B, with h = ixB. 
The behaviour of the energy for the cases considered above is sketched in Fig. 1. It is 
clear that the increase in magnetisation, when the temperature decreases, is slower in the 
case of the canonical T-distribution than for the conventional ensemble. This is because the 
latitudinal circles closer to the equator have larger weights than those closer to the poles. 
In Fig. 2 we plot the relationship between the specific heat and the temperature. 

We note that in the case of the canonical T-ensemble the heat capacity for this system is 
nonvanishing at zero temperature. For many bulk substances, on the other hand, the heat 
capacity vanishes as zero temperature is approached. Therefore, it would be interesting to 
enquire if a single electron has a different behaviour. To that end, we note that in the case 
of a general two state system, with energy eigenvalues Eq and Ei, the density matrix of the 
canonical T-distribution, in the basis of energy eigenstates, is 



+ 



+ 



(3{Eo-Ei) 1 - e-/3(^i-^o) 
1 1 



pap 



_(3{E,-Eo) ' l-e-^(^o-E,)_P''Pp- (22) 

This expression can be contrasted with the conventional result for such a system, given in 
formula (^Oj). In fact it is not difficult to obtain analogous expressions in higher dimensions, 
some of which may characterise other systems that can be implemented in experiments to de- 
termine whether the T-distribution introduced here provides a better account of phenomena 
in certain situations. 
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V. DISCUSSION 



The foregoing constructions are based on the quantum microcanonical postulate, which 
for an isolated system in equilibrium implies a uniform distribution over any energy surface 
of the quantum phase space. We therefore enquire whether this postulate can be derived 
from the basic principles of quantum mechanics. This would follow, e.g., if quantum dy- 
namics were ergodic on the energy surfaces £e- In the example of the spin one-half particle, 
ergodicity is indeed guaranteed by the periodicity of the Schrodinger evolution and the 
dimensionality of the energy surface. 

In higher dimensions, the generic energy surface £e is parameterised by (n — 1) angular 
parameters and n phase variables. For each fixed set of angular variables we obtain an ra-torus 
T"' in and by varying the angular parameters the energy surface is foliated by an (n — 1)- 
dimensional family of such tori. If we assume that the Hamiltonian is nondegenerate, and 
that the ratios Ei/ Ej {i ^ j) of the energy eigenvalues are irrational, then the Schrodinger 
evolution is nonperiodic on each T". As a result, ergodicity is guaranteed on the toroidal 
energy subsurfaces. On the other hand, unitary evolution does not change the angular 
parameters on Se- It follows that the energy surfaces are not fully ergodic with respect to the 
Schrodinger equation, and that we cannot deduce the microcanonical postulate directly from 
the principles of quantum mechanics. It would be interesting to see if this postulate could 
be derived as a consequence of a suitable generalisation of ordinary quantum mechanics, 
e.g., nonlinear or stochastic quantum dynamics. 

To summarise, we have formulated quantum mechanical analogues of the microcanonical 
and canonical ensembles by consideration of the geometry of the quantum phase space. In 
particular, a key distinction in our proposal is that phase space weightings are fully taken into 
account in the construction of equilibrium ensembles, whereas in the conventional approach 
such weights arise only in the event of energy eigenvalue degeneracies. We are thus able to 
give a coherent meaning to the quantum microcanonical ensemble, and derive an expression 
for the corresponding density matrix, a construction that is missing in the conventional 
theory. It may be that for bulk substances there is a further component to the dynamics, 
perhaps due to internal interactions at low temperature, that forces the quantum phase 
space to contract down to the lattice of energy eigenstates, with appropriate multiplicities. 
In that case we recover the conventional theory. It would be interesting to determine, by 
suitable measurements on quantum mechanical systems, whether the ensembles introduced 
here are indeed appropriate. 
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